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DEFINITION OF SYMBOLS 


Symbol 

A 


Anm 

dx 

d/dx , (') 
Ela 


EI(x) 

k 

L 


Lt 


rn 

rff 

n 

Un(x) 



^nO 
^■n 1 
^n2 

ir 


P 

<J>n 


Definition 

Peak stiffness variation, A = (Elmax - Elmin) / 2 
Factor defined by equation (A-2) 

Differential element length 

Denotes differentiation with respect to position 

Geometric average stiffness of beam design, Ela = (Elmax + Elmin) 12 

Bending stiffness function 

Dummy variable referring to mode number 

Span length 

Period of the stiffness function 

Dummy variable referring to mode number 

Mass per unit length 

Vibration mode number 

Denotes nth mode eigenfunction, bending 

Distance along span measured relative to support 

Dimensionless horizontal coordinate, x — x/L 

Kronecker delta function, equation (10) 

Stiffnesss perturbation parameters, e = A/EIa 
Dimensionless eigenvalue, \ 2 n = (mL 4 /EIa) co 2 n 
Zeroth-order eigenvalue solution 
First-order eigenvalue solution 
Second-order eigenvalue solution 
Pi 

Stiffness parameter, p = 1,2,3, i.e., the number of half periods p = 2 L/Lt 
Dimensionless Eigenfunction, 4>„ = Un/L 


v 



Zeroth-order eigenfunction solution 
First-order eigenfunction solution 
Second-order eigenfunction solution 
Denotes nth mode frequency, bending 
Partial derivative notation 
Integral sign 


Summation sign 



TECHNICAL REPORT 


DYNAMIC CHARACTERISTICS OF A VIBRATING BEAM WITH 
PERIODIC VARIATION IN BENDING STIFFNESS 


INTRODUCTION 


Vibrating beam theory has been considered extensively in the literature for any number of variable 
property states, ranging from structures with changing cross-sectional geometry to those of a composite 
nature. Solutions are obtained either in closed-form for a few simple cases, or they are pursued using 
numerical techniques. In the present study, a perturbation expansion technique applicable to continuous 
systems is used to develop a closed-form solution to the problem of a vibrating beam with bending stiffness 
periodic in the spatial coordinate. Results are compared to a finite element solution and verified experi- 
mentally using forced vibration of a test span. To the knowledge of the author, this specific beam problem 
(static or dynamic) has not been addressed in the literature. 

Application of periodic stiffness is recognized in the field of vortex-induced motion of transmission 
power lines [1]. In recent years a conductor, known as twisted-paired, has been developed that uses a vari- 
able diameter design to provide a changing conductor profile into the wind. Twisted-paired conductors are 
constructed by twisting together two identical standard round conductors with 360 deg twists occurring at 
set intervals along the span. The periodic nature of the twist causes a periodic variation in bending stiffness. 
Variable profile diameter results in non-uniform shedding of vortices, and hence excitation frequencies, 
along the span. Multiple vortex frequencies act to minimize wind energy transfer and detune vibration 
response. In conductor systems, influence of bending stiffness effects becomes extremely important in the 
vicinity of supports. 

The purpose of this report is to characterize the dynamic bending behavior of beams with periodic 
stiffness variation. Also, the models developed will provide insight into the behavior of similar type sys- 
tems with changing property states. 


EQUATIONS OF MOTION AND BASIC ASSUMPTIONS 


Consider the problem of the transverse vibrations of a straight beam with periodic variation in bend- 
ing stiffness along its length. The beam is assumed to be simply supported and long compared to its cross- 
sectional dimensions, and dynamic shear distortions and rotary inertia are negligible. We will also make 
the usual simplifying assumptions that Hooke’s Law holds and plane sections remain plane. Figure 1 shows 
a free body sketch of a differential element of the vibrating beam. 

We will proceed from the well-known differential equation of motion for the normal mode response 
of an undamped beam 


£ 

dx 2 


EI(x) 


d 2 U n (x) 

dx 2 


- mw 2 n U n (x) = 0 


( 1 ) 



where EI(x) is the bending stiffness function, m is the mass per unit length, U n (x) is the nth normal mode 
displacement, and co n is the nth normal mode frequency. The bending stiffness function is given by 


EI(x) = Ela + A cos 



( 2 ) 


where Ela is the geometric average stiffness of the beam design, Ela = [Elmax + EImin]/2; and A is the 
peak stiffness variation, A = [Elmax - EImin]/2. The period of the stiffness function is Lt. Figure 2 plots 
the function. Notice, at x = 0, the maximum flexural stiffness occurs, and the periodicity of the function is 
an even multiple of the span length. This particular stiffness function is characteristic of twisted-paired 
systems [1]. 

An equation of motion that models an undamped, vibrating beam with periodic bending stiffness is 
determined by combining equations ( 1 ) and (2) . The result is a fourth-order differential equation with vari- 
able coefficients. A closed-form approximate solution to this boundary value problem is obtained by using 
a variation of .the Rayleigh-Schrodinger expansion [2,3]. The solution in closed-form is extremely useful, 
since it clearly displays the influence of system parameters on response. Nayfeh [4] presents an application 
of a similar perturbation formulation for a simple linear second-order eigenvalue problem. 
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SUPPORT SUPPORT 



A dimensionless form of the governing motion equation is 


dx 


\V2 


1 + € COS p 7T xj 4> n " | - X 2 n c}) n = 0 


(See nomenclature section for a definition of terms.) 


and the corresponding boundary conditions for the case of simple supports are 


(3) 


4>n(0) = 4> n (D = o (4) 

and 

4>n"(0) = 4> n "(l) = 0 (5) 


The quantity [1 + e cos pirx] is the dimensionless bending stiffness. For convenience, the tildes are 
dropped in the remaining analysis. The coefficient e is a measure of the magnitude of the stiffness variation, 
e = A/EIa, and the parameter p is equal to the number of half periods of the stiffness function in a given 
span, p = 2L/Lt. 
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PERTURBATION EXPANSION SOLUTION 


The solution (4>,X.) of equation (3) is a function of the independent variable x and the parameters e 
and p. If the parameter e is equal to zero, the equation reduces to the case of a vibrating beam with uniform 
flexural stiffness whose eigenfunctions and eigenvalues are given, respectively, by 


4> n = V2sinmrx 

(6) 

X 2 n = n 4 rr 4 

(7) 


The V 2 " coefficient is arbitrary. It is picked so that the eigenfunctions are normalized according to the 
integral function, 


1 


J <j) 2 n dx = 1 
0 


( 8 ) 


The above eigenfunctions are orthonormal; i.e., 


1 

/ 4>n00 dx = 5 mn (9) 

0 


where, 8 mn , the Kronecker delta function is specified as 


8 


mn 


0 m 4 - n 

1 m = n 


( 10 ) 


When e does not equal zero, equations (6) and (7) are no longer valid and corrections must be added to them. 
An approximate solution is obtained by expanding both the eigenfunction and the square of the eigenvalue in the 
form of a power series in e; i.e., 


4>n — 4>nO + e4>nl + tj> n2 + ... 

(ID 

k 2 n — k n0 + e\„| + e 2 X n2 + ... 

(12) 



where 4> n0 and X n0 are the eigenfunction and the square of the eigenvalue when e equals zero; equations (6) and (7). 
An asymptotic expansion is generally valid only if e is small. By definition, the parameter e for beams with periodic 
stiffness may not be small, but it is always less than one. Corrections of the higher order terms are therefore neg- 
ligible, and the series converges eventually to the correct solution. 

Substituting equations (11) and (12) into equation (3) and equating like powers of e through order e 2 , we get 
the following system of differential equations: 


0(e°):<{> nO ,v -X n o<i> n o = 0 (13) 

0(e 1 ) : 4>n 1 1 v - X n0 4>n 1 = X n , 4> n0 - (d 2 /dx 2 ) [4> n0 " cos pirx] ( 1 4) 

0(e 2 ) : <j)n2 1 V - X^nO 4)n2 = ^n2 4>n0 + 1 1 - (d 2 /dx 2 ) [<j> n , " COS pTTX] (15) 


The problem of finding an approximate solution to equation (3) is now simplified to one of obtaining sequential 
solutions to equations (13), (14), and (15). To illustrate the procedure, the first-order correction is formulated in 
Appendix A. Equations (A-6) and (A-8) define the correction terms of the eigenvalues and eigenfunctions. 

Using equation (A-8) and recalling that 4> n = <f> n0 + e cf> nl where p = 1 ,2,3,..., the general eigenfunction 
solution 4» n is expressed in terms of dimensionless parameters for the case when n ^ p/2 as 



If n = p/2, the last term in the series is secular and hence vanishes. 

Using equations (7), (12), and (A-6), the eigenvalue solution X n is given in terms of the first-order correc- 
tion. The result is expressed in terms of the dimensionless variables as 

— n 2 TT 2 [l - (e/2)] 1/2 (17) 


Equation (17) is valid only for the vibration mode n = p/2. For all other vibration states the first-order correction 
term is zero. The second-order perturbation solution of the eigenvalue is determined using the same techniques 
previously developed for finding the first-order terms. Details of the formulation are outlined in Reference 1 . The 
general eigenvalue solution X n of the second-order expansion is given as 


— „ 2_2 

- n it 


1 - (e/2) - (e 2 /4) 


(n + p) 4 1 f (n-p ) 4 
_(n + p) 4 - n 4 J ^ 1 ^ .(n-p) 4 - n 4 _ 


1/2 


(18) 
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The first-order correction term in equation (18) is equal to zero for the vibration modes where n ^ p/2. For the case 
defined by n = p/2, the last term in the second-order correction is specified to vanish (i.e., this term is secular from 
4> n , solution). Notice that the eigenvalue \ n simplifies to the case of a beam with uniform stiffness when the perturb- 
ation parameter e is equal to zero. 

The general behavior of a beam with periodic bending stiffness variation is given by equations ( 1 6) and (18). 
As the vibration state approaches the anomaly occurring at n = p/2, the eigenfunction solution deviates from a 
simple sine wave displacement curve to a mode shape comprising other harmonics. For the case when n = p/2, the 
eigenfunction returns to the sine wave shape. The eigenvalue solution responds in a similar nature. At n = p/2 a 
jump in the eigenvalue occurs, since for this mode the harmonics of the stiffness function are secular. Results of the 
closed-form perturbation analysis have been checked using finite element results and results of experimentation. 


EIGENFUNCTIONS AND EIGENVALUES 


Effects of the perturbation terms on the eigenfunction solution are exemplified in Figure 3. The magnitudes 
of the bracketed terms in equation (16) are plotted as a function of the p/n ratio. Two distinct ranges are apparent; p/n 
< 1 and near p/n = 2. In these ranges, the perturbation effects are the strongest. At p/n = 2, a vibration state is 
defined where the lengths of the vibration loops match the period of the stiffness function. Based on the stiffness 
definition given, this mode defines maximum stiffness at the nodes of the vibration loops and minimum stiffness at 
the antinodes. At p/n = 1, both the maximum and minimum stiffnesses occur at the node positions in alternating 
sequence along the span. 

Figures 4 and 5 give eigenfunction comparisons of vibration modes in the general modal solution range. The 
vibration displacement amplitudes are normalized and plotted versus the normalized horizontal span coordinate 
(x/L). Recall, boundary conditions are simple support, and beam orientation at the supports is for maximum stiff- 
ness. Figure 4 plots eigenfunction solutions near the anomaly p/n = 2, where p = 64 and e = 0.4. The effect of 
periodic stiffness is to modulate the displacements of those vibration modes approaching the anomaly at n = p/2, or 
for this case mode 32. Similar displacement curves as patterned for modes 31 and 33 are characteristic of all vibra- 
tion modes near the anomaly. Close examination of the eigenfunctions reveals that the node (or antinode) locations 
are adjusting themselves along the span, and the longer vibrating loops result in lower midloop displacement ampli- 
tudes. Apparently, the beam attempts to minimize the elastic strain energy stored within the dynamic span by adjust- 
ing the lengths of the vibrating loops until the same average bending stiffness exists across each individual loop. 
Equalizing the loop stiffnesses may require the loops to have different lengths depending on the vibration mode, and 
a longer loop has greater mass. An equal partitioning of potential energy and thus kinetic energy between each of the 
loops results in lower vibration amplitudes for the longer vibrating loops. Loop stiffness calculations verify this 
reasoning. 

At the anomaly, modulation in the mode shape disappears, since for this case the individual loop stiffnesses 
are equal (i.e. , the lengths of the vibration loops match the period of the stiffness function). The same basic reason- 
ing holds true for the case where p/n = 1 (Fig. 5). Here also, the mode shape is sinusoidal - no modulation; and 
average loop stiffnesses are equivalent with maximum and minimum values defining the nodes of each loop. Addi- 
tional cases identified in Figure 5 are for small and large values of p/n. As p/n approaches zero, perturbations in the 
mode shape increase. At the opposite end of the scale, where the p/n parameter goes to infinity, the displacement 
shape is sinusoidal. Figure 6 characterizes the eigenfunction solution as a function of the perturbation parameter, e. 
As the magnitude of Epsilon (i.e., the stiffness variation) increases, the modulation effects become more 
pronounced. 
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VIBRATION AMPLITUDE 
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alt alts =lt> 









EIGENVALUE 


Eigenvalue solutions of the second order perturbation expansion for different cases of p are plotted in Figure 
7 versus the respective vibration mode. At the anomaly, p/n = 2, a jump in the eigenvalue occurs. The intensity of 
the jump increases with p. Physically, the jump identifies a sharp change in frequency (or stiffness) between 
characteristic modes of vibration. Mathematically, the jump is equivalent to removing the secular nature of the 
stiffness function from the eigenfunction solution. 

The closed-form perturbation expansion solution has been verified through comparisons with a finite 
element solution [1]. Although the findings are not formally documented herein, agreement between the analytical 
results is excellent. Some discrepancy does occur in the vicinity of the stiffness anomaly. This is apparently due to a 
sudden change in bending stiffness. Nevertheless, the qualitative picture remains the same. Two characteristic 
effects of periodic bending stiffness on dynamic response are determined: (1) periodic stiffness forces an anomaly in 
the system which results in a jump in the natural frequency, and (2) periodic stiffness acts to modulate the modal 
displacements in distinctive ranges of p/n. A qualitative explanation of the modulation and its effects on beam 
response is given in terms of energy principles. 


EIGENVALUE SOLUTION (e = .6) 



Figure 7. Eigenvalue solution for different cases of p. 
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EXPERIMENTAL VERIFICATION 


A series of tests were designed to investigate the dynamic response of beam type systems which have a 
periodic variation in bending stiffness. A stiff-string structure, known as a twisted-paired conductor, was the test 
candidate in the program. Table 1 summarizes the experimental test parameters and Figure 8 shows a photograph of 
the test span. Periodic variation in diameter profile of the twisted-paired conductor is compared to the uniform 
diameter of a standard conductor design. In stiff-string systems, elastic strain energy is stored in tension and bend- 
ing. If tension is constant along the span, then tension has minimal effect toward equalizing the variable flexural 
stiffness of the vibrating loops [1]. In other words, tension effects do not mask the effects of stiffness variation. 

Experimental data are compared to the finite element results for free vibration since the fixed boundary 
conditions are applicable. This type of boundary support keeps the end losses to a minimum. Internal damping of the 
conductor was also reduced by applying a high tension line force. The procedure of minimizing conductor system 
damping is necessary; higher harmonics are difficult to excite if mechanical damping is significant. The testing 
program used forced vibration response to study free vibration. If the span is tuned properly to a single natural 
frequency, contributions from all other harmonics are minimal. The vibration exciter unit was positioned near the 
span center to eliminate even harmonics from the general response. The added mass of the moving shaker element 
and span attachment fixture resulted in a shortening of the drive loop, and thus a lowering of its vibration amplitude. 
No attempt was made to decouple shaker mass from the conductor span. A V-scope attached to the center of each 
vibrating loop was used to measure the midloop amplitude displacements. The device is inexpensive and its 
accuracy is remarkable at 0.01 in. 



TABLE 1. SUMMARY OF TEST PARAMETERS 


Span Tension, T 9919 lb 

Span Length, L 144 ft 

Period of Stiffness, Lx 4.5 ft 

Stiffness Parameter 64 

Bending Stiffness, lb-in. 2 

El min. 386785 

El max. 2185202 

Ela 1285993.5 

Stiffness Parameter 0-7 

Weight, lbs/ft 1.164 

Fixed Boundary — 

Support Orientation Stiffness El min. 

Drive Unit MB C-Uc Shaker 

Weight of Moving Element of Shaker 2.20 lb 

Weight of Span Attachment Fixture 0.77 lb 



Typical test results are presented in Figures 9 and 10. Finite element displacement amplitudes are nor- 
malized using a method previously outlined [equation (8)]. Experimental amplitudes are normalized to one of the 
measured values - chosen in arbitrary fashion. Span length is used to nondimensionalize the horizontal coordinate. 
Since mode shape is symmetric about the span center, data results are shown only for half the span. Figure 9 gives 
the comparison for mode 27, an eigenfunction near the anomaly in the system occurring at n = 32. Although the 
comparison is not exact, the modulation in the eigenfunction response is proved physically to exist. The same 
general results are reported for all other modes near the anomaly. Figure 10 compares the experimental and analyti- 
cal data of mode 19, an eigenfunction well removed from the stiffness anomaly. Agreement between the analytical 
and test data is excellent and the mode shape is sinusoidal. Some discrepancy does occur in the node positions near 
the drive location where the measured loop lengths are shorter. Shacker attachment changes the stiffness and mass 
of the drive loop. Shifting of the nodes tends to compensate for these effects. 


CONCLUSIONS 


The dynamic response of beams with periodic stiffness contrasts significantly with the vibration behavior of 
standard beams. Linear vibration theory was used to develop a stiffness model and characterize response. Using a 
perturbation expansion, a closed-form solution of free vibration was formulated for the case of simple supports and 
periodic stiffness variation. The technique worked exceptionally well when the stiffness parameter was slowly 
varying. Applications of variable tension, mass, and area are natural extensions of the theory. The main conclusions 
are summarized below. 

1 . Periodic bending stiffness forces an anomaly in the system which corresponds to the vibration state 
where the loop length matches the period of the stiffness function. Physically, the anomaly denotes the vibration 
mode for which loop stiffness changes most rapidly. The result is a jump in natural frequency. The perturbation 
solution loses some accuracy for those vibration modes near the anomaly; however, the qualitative characteristics of 
the response remain the same. 

2. The stiffness parameter acts to modulate the modal displacements in two distinct ranges of vibration: p/n 
< 1 and near p/n = 2. Experimental evidence is presented which supports these findings. 

3. Modulation in modal displacements is explained in terms of energy principles. The beam attempts to 
minimize the elastic strain energy stored within a dynamic span by adjusting the lengths of the vibrating loops until 
the same average bending stiffness exists across each individual loop. Equalizing the loop stiffnesses may require 
the loops to have different lengths depending on the vibration mode, and a longer loop has greater mass. An equal 
partitioning of potential energy and thus kinetic energy between each of the loops results in lower vibration ampli- 
tudes for the longer vibrating loops. 


12 
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DISTANCE ALONG 1/2 SPAN (X/L) 

Figure 9. Comparison of eigensolution with experimental test valves, mode 27. 


MODE 19, C= .7, p = 64 



0 .25 5 

DISTANCE ALONG 1/2 SPAN (X/L) 

Figure 10. Comparison of eigensolution with experimental test valves, mode 19. 





APPENDIX A 


FIRST-ORDER CORRECTION TO THE PERTURBATION SOLUTION 


The first step in the first-order correction to the perturbation solution is to substitute the zeroth-order solu- 
tion, equations (6) and (7), into equation (14). After taking the appropriate derivatives and using trigometric identi- 
ties, the result simplifies to 


4> n , 1 v - n 4 tt 4 4> n | = V2 {X n | sin nrrx - (n 2 ir 4 /2) (n + p) 2 sin [(n + p)rrx] - (n 2 iT 4 /2) (n-p) 2 sin [(n-p)7rx]} ( A- 1 ) 


Next, assume that the solution 4> n i can be expressed as a linear combination of the zeroth-order eigenfunctions <j) n0 ; 


CO 

<t>nl = V2S A nm 
m = 1 


sin rmTX 


(A-2) 


This solution satisfies the boundary conditions on 4> nl . Taking derivatives and substituting into equation (A-l) 
yields 


00 

2 -n- 4 (m 4 - n 4 ) A nm sin mirx = \ n | sin nirx - (n 2 Ti 4 /2) (n + p ) 2 sin(n + p)itx 
m = 1 

- (n 2 orr 4 /2) (n-p) 2 sin(n-p)TTX (A-3) 


Multiplying equation (A-3) by sin k-nx and integrating from 0 to 1 using the orthonormal property, equation (9), we 
obtain 


7r 4 [k 4 -n 4 ] A nk = \ n) 8 nk - n 2 ir 4 (n + p) 2 J sin(n + pkrrxsinkTTxdx 

- n 2 tt 4 (n-p) 2 f sintn-p^xsinkTTxdx (A-4) 


If k = n, the left-hand side of equation (A-4) vanishes, hence 


^-n I 


1 

n 2 tt 4 (n + p) 2 / sin(n + p)rrx sin n TTxdx 
0 

1 

+ n 2 7T 4 (n-p) 2 / sin(n-p)TTXsin n ttx dx 
0 


(A-5) 
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The above integral expressions evaluate only when n = p/2, that is when the vibration mode number corresponds to 
the span length-stiffness function ratio (n = L/Lt). This is the anomaly that makes the periodic stiffness problem so 
interesting. 

Equation (A-5) then calculates the eigenvalues of the first-order expansion as 


\ n , = -n 4 7r 4 /2 when n = p/2 


(A-6) 


Note, \ n) = 0 for all other values of n. The above condition removes the secular terms from the solution when k = 
n. If k ^ n, then equation (A-4) simplifies to 


1 

Ank [-n 2 (n + p) 2 /(k 4 -n ')]/ sin (n + p) T7X sin krrxdx 

0 

1 

+ [-n 2 (n-p) 2 /(k 4 -n 4 )] J" sin (n-p) ttx sin k-nx dx (A-7) 

0 


Because of the stipulations on the parameter p, Ank in equation (A-7) calculates non-zero values only for the two 
cases: (1) k = n + p and (2) k = n - p, where n £ p/2. The general solution of the eigenfunction 4> nl , given by 
equation (A-2), is finally expressed as 


-n 2 (n + p) 2 
= VT [(n + p) 4 -n 4 ] 


sin (n + p)rrx 


. n 2 (n— p ) 2 

V2 [(n-p) 4 -n 4 ] 


sin(n-p)-JTX + V^AnnSinmrx 


(A-8) 


For the vibration mode corresponding to n = p/2, the second term in equation (A-8) vanishes, since for this mode it 
is secular in nature. Keeping this in mind we can say that the cj> nl solution is valid for all vibration states where p = 
1,2,3,... . The coefficient A nn is determined by the normalizing function 

1 

/2 4. nO <f> nl dx = 0 (A-9) 

0 

For this case, A nn is calculated to equal zero. 
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